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Approximating Linearizations 
for Nonlinear Systems 

L.R. Hunt, R. Su, and G. Meyer 
ABSTRACT 

Given a nonlinear control system 


m 

x(t) = f ( x ( t ) ) + l u . (t) g . (x (t) ) 

i=l 1 1 


on E n and a point x^ in !R n , we want to approximate the system near x^ 
by a linear system. Of course, one approach is to use the usual Taylor 
series linearization. However, the controllability properties of both 
the nonlinear and linear systems depend on certain Lie brackets of the 
vector field under consideration. This suggests that we should con- 
struct a linear approximation based on Lie bracket matching at Xq. In 

general, the linearizations based on the Taylor method and the Lie , 

.. . > 

bracket approach are different. However, under certain mild assump- 
tions, we show that there is a coordinate system for near in 
which these two types of linearizations agree. We indicate the im- 

f 

J 

portance of this agreement by examining the time responses of the ndn- 
linear system and its linear approximation and comparing the lower * 
order kernels in Volterra expansions of each. 
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I. Introduction 

Suppose we have a nonlinear control system that we wish to sim- 
plify in some way. An approach that has received much attention in 
the literature is the exact linearization, whereby the nonlinear sys- 
tem is transformed to a linear system. Both theoretical problems [1] , 
[2] , [3] , [4] , [5] , [6] , [7] and practical applications [8] , [9] , [10] , [11] , 
[12] , [13] , [14] , [15] , [16] , [17] concerning this method exist. Theoret- 
ically, a nonlinear control system in which the controls enter linearly, 
is "(locally) equivalent" to a controllable linear system if and only 
if (i) a certain set of vector fields is linearly independent and (ii) 
related sets of vector fields are involutive (see [3] and [5]. These 
two conditions will be explicitly stated later. 

The ideas in this paper were generated by the desire to construct 
an approximate transformation (if an exact one cannot be found) for a 
nonlinear system that is transformable to a linear one. Meyer [12] 
used linear Taylor series expansions about certain points along a 
flight trajectory to discover approximate transformations. However, 
this type of linearization did not in general reflect the rich differ- 
ential geometry inherent in the assumptions (i) and (ii) previously 
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mentioned. Therefore, we introduced in [20] a linear approximation 
based on the important Lie brackets. 
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Given a nonlinear system which satisfies condition (i) (and not 
necessarily condition ( i i ) ) and a point x Q in state space, there are 
two types of linearizations of interest. One depends on Taylor series, 
the other on Lie brackets, and in most cases these do not agree. How- 
ever, if Xq is an equilibrium point of the drift term of the system or 
if we are operating about a known trajectory containing the point x^, 
then the two linearizations are the same. 

Recently, the authors have discovered for nonlinear systems sat- 
isfying the condition (i) that there exists a coordinate system (called 
the s coorinates) on state space in which the nonlinear system is in 
a particularly nice form [18], [19]. This state space coordinate 
system appears in the literature in [1] and [3], but the canonical 
system expansion of a general nonlinear system presented in [19] is 
new. 

The main point of this paper is to show that the two kinds of 
linearizations coincide in the s coordinate system that we studied 
in [19]. We consider formal Volterra series and stress the importance 
of this agreement in computing Volterra kernels. From an input to 
output (input to state in our case) point of view, for any appropriate 
input the error between the time responses of the actual system and 
the approximating system (in the special coordinates) propagates 
like 0 ( | t | ^ ) in the single input c&se. 

We show how to compute the coordinate changes to move from the 
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original coordinate system to the one in which the linearizations 
agree, but it is not always possible to do this in practice. However, 
we are still able to find the linear part about of the nonlinear 
system in the special coordinates. 

In Section 2 of this paper we present definitions, review the 
exact linearization results, and indicate the desired coordinate 
system for a nonlinear control system. Section 3 consists of an 
example, the main result, and the interpretation of this result in 
terms of Volterra kernels. 


II. Preliminaries 

For C°° vector fields f (x) and g (x) on B n we denote the Lie 
bracket (this is the negative of the usual definition) 

_ 9g f 

by [f • where and are Jacobian matrices. We also define 

(ad°f,g) = g 
(ad 1 f,g) = [f,gj 

(ad k f ,g) = If , (ad k-1 f ,g) ] . 

Let h(x) be a C°° real-valued function on B n and 


f (x) = 


n 



4 


a C°° vector field. The Lie derivatives of h with respect to f ,L^h(x) , 
is defined to be . 


<dh, f > 



+ 


ah f 

3x 2 2 


+ . 


. + 


3h f 

8x n* 
n 


We take 


L^h(x) = h(x) 

L^h(x) = L^h(x) 

L^h(x) = <dL^ 1 h(x) / f>. 


Given a nonlinear control system on H n 

m 

(1) x(t) = f(x(t)) + l u. (t)g. (x(t) ) , 

i=l 1 1 

where f,g, , ...,g are C°° vector fields, and a point x n in F n , we are 
1 in u 

interested in finding linear (affine) approximations of the form 


(2) 


x ( t ) = f 


<V ' 


Ax_ + Ax + 


m 

l 

i=l 


u . (t) b . 
x i 


which are appropriate for use in control problems. 


Let • k 2 ' * * * ' K m positive integers such that ^2— * * *— K m-^ 

and k 1 +ic = n. We take the sets 


Kf-i 

c = (g^ , ff ,g 1 ] , • . . , (ad f ,g.) ,g 2 , [f ,g 2 J » • • • » 

< 2 “! <^-1 
(ad f ,g 2 ) , . . . ,g m , [f ,g m J , . . . , (ad m f,g m )l 
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k .-2 


j ^ £ f Q r *3 2 f ^ ^ ^ ‘ m * ' 

K . K 

(ad ^~ 2 f ,g 2 ] , . . . ,g m , [f ,g m ) , . . . , (ad ^~ 2 f,g m )} 


for j=i , 2 , . . . ,m. Suppose near system (1) satisfies the two condi- 
tions : 


(i) 

(ii) 


The set C spans n dimensional space and the span of C. 


equals the span of C^fl C for each j=l,2,...,m. 

Each set , j=l , 2 , . . . ,m, is involutive; i.e. the Lie bracket 
of any two vector fields in Cj is a linear combination of 
the vector fields in C . • 


D 

Then it is proved in [3] and[5] if f(Xg)=0 and [21] if f(Xg)^0 that system 
(1) is locally equivalent to a controllable linear system with Kronecker 

indices K 2. ' *2 ' * * * ' K m ^ we can renum ^ er t ^ ie to ma ^ e 

>^<2 > • ♦ K m -if necessary). Hence there are new state and control 
variables in which system (1) is actually a linear system for x near 
Xq. This is called an exact linearization of (1) . 

If the above two conditions hold and the state and control trans- 
formations can be found (a method for constructing such transformations 
is described in 15]), then it is not necessary to approximate the non- 
linear system (1) by a linear system because in the correct coordinates 
it is a linear system. 

Suppose we assume that assumption (i) holds for system (1) but 
discard condition (ii) . We present a coordinate system (called the 
s coordinates) on B n near x^ in which our nonlinear system takes a 
particularly nice form. In fact, if condition (ii) also holds 
we have a pure feedback system as in [19] (related results are in 
[3]), and it is trivial to move from this form to a linear system. 



In a special case these pure feedback systems are called block 
triangular by Meyer and Cicolani [ 9 . 1 . 
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Emphasizing that we are working under condition (i) only, we re- 
order the elements in the set C to reflect descending orders in the 
superscripts on the ad's and ascending orders in the subscripts of 
the g's. We call this reordering C' and the first element of C' is 


< 1-1 


< 2 -l 


(ad f,g^). If k ^ =k 2' t * ie secon< ^ element of C* is (ad f,g 2 ), 


<^-2 

and if <^ ><2* it is (ad f,g^). If <-^=< 2 = K 3 ' t * ie third element 

<2~1 c^-2 

is (ad f ,g^) / if <^= <2 ><3, it is (ad f,g^) f and if 

< 1~3 < 2 ~1 

it is (ad f,g^) or (ad f,g 2 ) depending on whether k^-1>k 2 or n °t 

The process can be continued, and the last element in C' becomes g . 

^m 

This order is simply the opposite lexicographic order. 

For x Q =(x 10 ,x 20 , . . . ,x n() ) given and s Q = (s io' s 20 ' * * * ' S n0 ) we solve 
in order the following systems of ordinary differential equations with 
initial conditions; 


dx(s 1 ) < 1 ~1 

-j- (ad f, 9l ), x(s 10 ) = x ( 


( 3 ) 


dx (s. ,s„) 

^ = 2nd element of C',x(s^,s 2 q) = x(s^) 


dx(s. ,s_,s 3 ) 
cfiT 


= 3rd element of C' ,x (s^, s 2 , s^q) = x (s^ , s 2 ) 


dx ( s . , s 2 > • • • / s ) 

ds ~ = ^m' X ^ S l ,S 2 ' ' * ' ,S n-l' S nO^ = x < s i ' s 2 ' * * * ' s n -l ) 


n 


Since the set C 1 consists of linearly independent vector fields, by 
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the inverse function theorem we can solve for s.,s_,...,s as func- 

1 z n 

tions of x.,x_,...,x near x n . Moreover, we can take the point 
l i n 0 

Sq = 0, the origin in s-space. 

We now view the s coordinates geometrically and introduce a 
sequence of manifolds Sq/S^,..., 5 in the following manner. is 

the point 0 and is the one-dimensional integral manifold of 

Kf-l 

(ad f,g^) through the point 0. Similiarly, is the two- 
dimensional manifold constructed by taking the integral curves of the 
second element of C' through S^. Likewise, is the three-dimensional 
manifold formed by merging the integral curves of the third element of 
C' through S^. Continuing in this manner, the manifold S n is a neigh- 
borhood of the point 0 in 3* n on which the above process is guaranteed 
from the inverse function theorem. Hence in the s coordinates 

g is the vector field 


g , (or the second to last element of C') 
m— 1 


0 

0 

I when restricted to S . , 
. n— i 

0 

1 

0 . 


g (or the third to last element of C') 
m-z 


is the vector field 
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0 

0 


is the vector field 


I when restricted to S - , 

. n-2 


1 

0 

_0 


"0 

1 

0 


the second element in C' 


is 


when restricted to S 2 » and 


Loj 


K 

(ad 




on S i . 


0 

0 


The above notation assumes that we have written the nonlinear 
system (1) in the s coordinates as 


( 4 ) 


m 

s = f (s) + l u.g . (s) , 
i=l 1 1 


^ 9 X 

where the new f is (-r— ) f(x(s)) and each new g. is (y— ) g. (x(s). 

o S 1 o S X 

We remark that in the s coordinates. 


a 

^m 3 s 

n 

9 

g , = on S , 

^m-1 3s . n-1 

n-1 


K 


1 


-1 

f,g 1 ) 


3 


3s 


1 


1* 


(ad 


on S 
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These facts will prove to be useful in our later work. 

We shall return to the nonlinear system (1) and introduce two 
kinds of approximate linearizations for nonlinear systems. 

III. Approximate Linearizations 
Given system (1) and point x^, we suppose that condition (i) is 
satisfied. We can do the usual Taylor series approach to find the 
linear approximation (about Xq and zero controls) 

m 

(5) x(t) = f(x n ) - Ax + Ax + l u . ( t ) b . , 

u u i=1 i i 

8 f 

where A = (x.) and b. = g.(x n ), i=l,2,... f m. This is called the 

OX U 1 1 U 

tangent model in [12] . Using Lie bracket matching at x^ we arrive at 
the modified tangent model [20] 

m 

(6) x(t) = f(x n ) - Ax n + Ax + l u. (t)b. , 

u u i=l 1 1 


where A,b, ,...,b are defined by 
1 m 



A k 

A b^ = 

(ad k f ,g 1 ) (Xq ) , 

k 0/lf • • 

(7) 

.k 

A = 

♦ 

(ad k f ,g 2 ) (x q ) , 

k — 0 , 1 , . . . , <2 


k 1 
A b = 
m 

(ad k f,g m ) (x Q ) , 

m 


In general the A matrices defined by (5) and (7) are different. 

The advantages and disadvantages of each of these two types of 
1 ineari za t ions will be stressed in terms of the formal Volterra series 
introduced later. We now show that for classical problems in control 
theory, these linearizations agree. 
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* 

Suppose Xg is an equilibrium point of the system x=f (x) in (1) 

(i.e. f(xQ)=0) and assume Xq= 0. The system (5) given by Taylor series 
has the property that A b i =(ad f,g^)(0), k=0 , 1 , 2 , . . . , , and i=l,2,...,m. 
Thus the tangent model and modified tangent agree in this case. 

Suppose ip is a trajectory of system (1) corresponding to all u^=0; 
in other words <p=f(tp(t)). We let 

m 

z = A(t) z + l u . (t) b. ( t) , 

i=l 1 1 


where 


z = x-cp(t) 

A = HMtl) 

b i = g i (tp ( t) ) . 

Setting (p(tg)=XQ and r=A(t)- a time varying Lie derivative, we 

find (as Hermes did in [22]) that 


W 


r W 


r2b i (t o> 


W - i - 




lx* x o’ 9 i (x O ) 


tf ,g i J (x o ) ' 


1^2^ • • • f m / 

3b i (t 0 ) 

at 

a g . 

" -aT (x o )f(x o ) 

i 1,2,..., m , 


(ad f tg/) (Xq ) , i = 1, 2, 


,m, 


Hence linearizing about the trajectory 4 >(t) and evaluating at XQ=ip(tQ) 
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show that the Taylor series approach and Lie bracket matching method 
yield the same result. 

3 

We now present an example in P to show that these approximating 
linearizations can be different. 


Example 3.1, 


On P' 


we take the single input system 


-1 


- 


— -i 

*1 


X 2 +X 1 X 3 


X 3 


= 

X 3 

+ U 

0 

*3. 


X 

L 


- 1. 


f (x) +ug (x) 


and Xq= (x^q ,x^q ,x^q) • The tangent model in this case is 


• 


- 


- 


- 


- “ 



X 1 


X 20 +X 10 X 30 


x 30 

1 

X 10 


x l _x 10 


X 30 

X 2 

= 

X 30 


0 

0 

1 


X 2" X 20 

+ u 

0 

_ X 3_ 


X 

0 

1 


1 

0 

0 __ 


- X 3" X 30- 


_ 1 


= f(x Q ) + A (x-Xq) + ub. 


Computing we find that 



"2 1 


"2 


X 30 +X 10 


x 30 

Ab = 

1 ] 

, [f ,g] (X q ) = 

1 


„ X 30 


1 

X 

u> 

o 


and if x^^O, the tangent model and modified tangent model cannot 
coincide. 

The work of Hermes [ 22 ] on controlling a system along a trajectory 

indicates for this example that the important Lie brackets are 

2 

g(x Q ) , [f ,g] (x Q ) , (ad f ,g) (x Q ) . These are linearly independent at any 
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point Xg. Thus if we are operating near a point x^ with x^^O , the 
modified tangent model seems appropriate. On the other hand, if the 
control u=0, the standard linear Taylor series approach seems 
reasonable. In our later discussion of Volterra series, these 
observations will be explained. We remark that in our s coordinate 
system, we do not have to choose between these types of linearizations. 
The authors wish to thank a reviewer for shortening the proof of the 
following result. 


Theorem 3.1. Given the nonlinear control system (1) satisfying con- 
dition (i) near the point x^elR 11 , there exist a local coordinate 
system on !R n at Xq for which the tangent model and modified tangent 
model agree. 


Proof . The s coordinates are the obvious candidates so we assume 
our nonlinear system is given by equations (4) 

. m 

s = f (s) + Z u.g. (s) . 

i=l 11 

First we construct the tangent model at the point 0 where x(0)=Xg. 
We write 


f(s) — f (s. ,s 0 , . .. ,s ) 
12 n 


and expand in a Taylor series to find 


f ( s ) = f (0 ) + ? |f- (0) s + 0 (s 2 ) 

i=l 5 i 


Since g is equal to — , g . is equal to -r on S . 

^m 9s m-1 9 s . n-1 

n n-1 


# • • * 9 


V 1 9 

(ad f,g) is equal to jg- 


on S we have 
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If (01 = (f,g ro )(0) 


3f 

3s 


n-1 


(0) - C f ,g m _ 1 J (0) 


3 f ^ l 

(0) = [f , (ad 1 f,g L ) ] (0) 


(ad "r ,g^) (0) 


Expanding 


m 


s=f(s)+ Z u.g.(s) 
i=l 1 1 


in a Taylor series with zero controls yields the tangent model 

m 

s = f (0) + As + Z u.b., 

i=l 1 1 

9 f 

with b. = g.(0),i = l,2,...,m and A = (0). Then equations (7) 

11 d S 

are easily verified for A and b.. □ 

A discussion of formal Volterra series, in which questions con- 
cerning convergence are ignored, is appropriate. We take the non- 
linear system (1) and add as output the identity function on I* n , 
i.e. 


y = h (x) = (h 1 (x) ,h 2 (x) , . . . ,h n (x) ) 
h^x) = x 1 
(8) h 2 (x) = x 2 


h (x) = x . 
n n 

If we are concerned with convergence, then we must take f / g«,...,g m 
to be real-analytic. However, since we are interested in low order 
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Volterra kernels, we stay with the C assumption and consider ex- 
pansions of low order plus remainder as in [23]. 

We assume for the rest of our discussion on Volterra series that 
the initial value problem x=f (x) ,x(0)=Xq has a solution on [0,T] , for 
some T >0, and the real-valued inputs considered are in L -([0,T]). 

This allows us to discuss finite Volterra expansions with remainder 
(analogously to that of Taylor series) as in [23] . The time t will be 
restricted to the set [0,T]. 

From [24] we take the Volterra expansions for the system (1) with 
output y=h(x) 


(9) 


y (t) 


w 0 (t) + 


m 


l 

i=l 


t 

w i (t, x^) u.^ (x^di 

0 



T 2 )u i 1 <T l ,U i 2 (t 2 )dT l dT 2 +0(t3) 


where 


(using <dh,f> = 


<dh 1 ,f> 
<dh 2 , f > 


<dh n ,f> 


= L^h, etc.) 


w 0 (t) 


I^M X ih = etfh lx 

k=0 1 0 K * X 0 


(10) w i (t# T 1 ) 


00 k 0 k 

1 L f L g L f h l X 

k 2 ,k 1 =0 g i 0 


k k 2 
(t - T l ) T 1 

-kT^y 


# i-1 $ 2 f • • • f in 9 



15 


w . . 

1 1 1 2 


(t, T 


1' T 2 j 



^1 ^2 k "? 

(t - T i> (T r T 2> t 2 

k l !k 2 !k 3 ! 


L l fi 2 


= 1 , 2 , 


,m. 


Here the notation from [24] has been extended to multi-input systems, 
and the infinite series are to be taken formally and will be truncated 
in our discussion. 

Now 

Lfh (Xq) = h (Xg) = Xq , 

(11) L^h (x Q ) = <dh,f>(x Q ) = f(x Q ) 

L f h <V - = M (x 0 )£( V- 


Recall that the important consideration in this paper is the approxi- 
mation of the nonlinear system (1) by a linear system of the form (2) . 
Suppose we consider Vol terra expansions in the form (9) for systems 
(1) and (2) , both having the identity as output. If (2) is formed by 
the usual Taylor approach (i.e. we have the tangent model) , then for 
the first three terms in w^(t) for the system (2) we obtain 


( 12 ) 


L f (Xq ) +A(x-x q ) h * x 0 ^ 

L f (Xq) +A (x-Xq) h(x 0 } 
2 

L f (x Q ) +A(x-x q ) h (x 0 } 


x 0' 


f {x 0 ) 


Af (x Q ) 


ff tx o )f(x o> 


Hence we have agreement in the kernels w^(t) for the nonlinear system 
and the tangent model through order t terms. This is a significant 

characteristic of the linear Taylor series expansion for our nonlinear 

3 f 

system. If A in (2) is not -^-(x^), as can occur in the modified 
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2 

tangent model, this agreement through order t is not assured. 

t f * 

Letting e h| denote the flow of the system x=f (x(t)) starting 

X 0 

at Xq, we can rewrite the kernels w^(t,i^) as (see [24]) 


(13) 


w.(t,T ) = l (ad k r f,g.)e tf h | , i = 

1 1 k=0 K ’ 11 x 0 


We emphasize the appearance of the Lie brackets (ad f,g^) in these 

kernels. The modified tangent model appears to be more natural than 

the tangent model because of this Lie bracket matching at x^ through 

order k=K . in w . . 

l i 

It should be obvious that by working in the s coordinate system, 
where the two types of linearizations agree, we have nice approximation 
from the input to state map point of view. Assume that system (1) is 
in the s coordinates (i.e. let x=s) and suppose y and y L are the 
Volterra expansions for systems (1) and (2) respectively. Then, since 
the Taylor approach and Lie bracket matching method agree in (2) , 


_ t 
in t °° 


y-y L = 0(t 3 ) ♦ I f l ±£l((adV, gi ,e« h l x 

i=l ! W 11 *0 1 1 


k.k._ t(f(x 0 )+A(x-x 0 ) 


(14) 


m r T l 

lx„ ,u i <T i )dT l + 1 l . w i.i, (t l T l' T 2 )u i, (, l >u i ‘ T 2 )dt l dT 2 +0< l t l 


Here the w. . (t.T.,T_) are for system (1), and the corresponding ker- 

^2 X i. 4 

nels for the linear system (2) are" of course zero. We are interested 

2 • • 

in those terms that contribute to degree t or less, the remaining 

3 1 

terms being moved to 0(|t| ). Hence we consider £ in (14), and in 

k=0 

fact examine only k ^=0 and k^-O, k^ = l and k2=0, k ^=0 and k2=l for 
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w. in (10) and k =k =k =0 in w. . 
i 123 1 1 1 2 

Computing we find 
t 


m f . 3g. 

y ~ y L = J. J (g i (x o ) + 5 (x o ,9 i (x o ut ' T i l+ -3F (x 0 )f(x 0 lT i' b i‘ Ab i (t ‘ T l ))u(r l l 

0 


_ t Tjg. 

m r r 1 


+ . I 

x l ' x 2 1 o 0 


T3T (x 0 )g i , (x 0 )u i, lT l )u i < T 2 )d, l dT 2 + 


(15) 


m t 

= -l 
i=1 0 


3f 


( (- If # g ± ] (X 0 ) + Ab i )T 1 + (U(X 0 ) gjL (X 0 )- Ab^ (t))u(T 1 )dT 1 


« rYl^i. 3 

. I -j J “13T (x 0 )g i. (x 0 )u i 1 (T l )u i, (T 2 )dT l dT 2 + 0( l t > U 


1 ' 2 0 0 


Using the fact that the tangent model and modified tangent model agree 
we obtain 


(16) 


y-y L - 


m r 1 r T 1 3g i, 

. 1. J J -33T (x O )g i < x 0 )u i 1 (T l ) % (T 2 )dT l dT 2 +O( l t l 

X 1 ,X 2~ X 00 l L t 


Suppose that we have a single input system (i.e. m=l) . In the s 
coordinates, g^ is the constant vector field 


and y-y =0 ( | t | ). For a two input system, 
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Hence 

(V 1 "l 

y -y L = l j n -(x o ) (gi( x 0 )u i (T 1 )u i (T 2 )+g 2 (x 0 )u l (T l )u 2 (T 2 ) ) dT i dT 2 + °( l fc 
0 0 


However, if and g2 are both constant vector fields in the s coordinate: 
(e.g. this can be done if then y-y L =0(|t| ). 

Thus a pattern emerges which can be extended to a system having 
any number of inputs. The importance of the s coordinate system (and 
of the agreement of the tangent model and modified tangent model in 
these coordinates) in time response studies has been proved. 

Starting with system ( 1 ) satisfying condition (i) in any x coordi- 
nate system, how do we find the tangent model, and thus the modified 
tangent model, in the s coordinate system? It certainly is not always 
possible to solve in closed form the systems of ordinary differential 
equations (3) . 

m 

Given x = f (x) + £ u.g(x) and a point x n =x( 0 ) we have 

i=l 1 u 


s = f(x(s)) + l g (x (s) ) . 


From ( 3 ) we know tt— where the entries are functions of x. ,x_,...,x . 

dS 1 1 2 n 

'N v <"\ " i 

Since -5— is invertible we obtain ( 4 ^-) with entries as functions of 

d S d S 

X l' x 2 ' * * * ' x n * T ^ e tan 9 ent niodel at 0 is 
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8 x 


-1 


5 ~ ( 8? } (x 0 )f(x 0 > 


8 x 


-1 


31(53) f(x))(x 0 ) 

+ 53; < 3 ?> <* 0 > s 


m ~ -1 

+ J 1 u i (< Ti ) (x o )g(x o )> * 


The paper [20] is written from the point of view that the modi- 
fied tangent model is more natural than the tangent model for con- 
structing approximate transformations to linear systems for exactly 
linearizable nonlinear systems. However, since these two models 
agree in the s coordinates, no choice need be made. We simply find 
the tangent model in the s coordinates and apply the approximate 
transformation theory of [12]. In designing a trajectory autopilot 
for VSTOL aircraft, the method of [12] has been successfully tested 
in flight simulation. 

In this article we have considered two types of linearizations 
of a nonlinear system about a point x^. We have found a coordinate 
system in which these agree and have shown the value of this in 
examining input to state time response through Volterra expansions. 
Some of the results of this paper are presented in preliminary form 
in [25]. 

Recent results by Krener [26] on approximate linearization by 
state feedback and coordinate changes are quite interesting. 
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